Distribution of Eigenvalues for the 
Modular Group 

E. Bogomolny * F. Leyvraz "I" C. Schmit 

Division de Physique Theorique-'- 
Institut de Physique Nucleaire 

91406 Orsay Cedex, France 



Abstract 

The two-point correlation functions of energy levels for free motion on 
the modular domain, both with periodic and Dirichlet boundary con- 
ditions, are explicitly computed using a generalization of the Hardy- 
Littlewood method. It is shown that in the limit of small separations 
they show an uncorrelated behaviour and agree with the Poisson distri- 
bution but they have prominent number-theoretical oscillations at larger 
scale. The results agree well with numerical simulations. 
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1: Introduction 



Free motion on constant negative curvature surfaces (CNCS) generated by dis- 
crete groups is the oldest and in some sense tlie best example of classically 
chaotic motion (see e.g. |l| ^). In recent years this subject has attracted 
wide attention also within the context of quantum chaos. There the main ques- 
tion is the way in which classical chaos manifests itself in the properties of the 
corresponding quantum systems (see e.g. |^). An important property of 
CNCS models is the existence of an exact relationship between the density of 
eigenvalues of the Laplace-Beltrami operator on the surface (= energy levels) 
and the geodesic on the surface, which correspond to classical periodic orbits. 
This is known as the Selberg trace formula (see e.g. (6[ 0). For arbitrary 
systems, only an approximate connection of this type is known, namely the 
so-called Gutzwiller trace formula ^, which is asymptotically valid in the 
limit of highly excited states but does not have a good estimate on the error. It 
is therefore the coexistence of hard classical chaos and the exact Selberg trace 
formula which makes the study of CNCS models so important. 

The simplest hallmark of classical chaos for a quantum system is the na- 
ture of the spectral fluctuations of the energy levels. It was conjectured that, 
for ergodic systems with strong chaotic properties, the fluctuation properties 
of energy levels should be that of the classical random matrix ensembles 
This result has found considerable numerical confirmation, but numerical work 
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on various models on CNCS jn], |l^, |l^, Q gave unexpected results. It was 
observed that the distribution of energy levels for these systems is quite close 
to a Poisson distribution, normally typical of integrable systems and not to 
random matrix ensembles which are characterized by strong level repulsion. 

In 17 1, it was shown that this anomalous behaviour could be traced 
back to a non-generic feature of these systems, namely the fact that the under- 
lying group belonged to a very specific subclass of the discrete subgroups of the 
motions of CNCS, namely the so-called arithmetic groups. 

Arithmetic groups are groups which permit a representation by 71 x 71 ma- 
trices with integer entries pO| . The important consequence of the arithmetic 
nature of these groups is that in such cases the corresponding CNCS shows an 
exponential proliferation of geodesies having exactly degenerate lengths |l^ . 
It is the cumulative effect of the interference of these degenerate orbits which 



leads to the Poisson-like distribution of energy levels. In 16, ^ the two-point 
correlation function of energy levels was computed in the diagonal approxima- 
tion , where one neglects all correlations between orbits that are not exactly 
degenerate. It was shown that this function definitely differs from the result 
predicted by random matrix theory pO[ |. Unfortunately, the diagonal assump- 
tion is quite crude and is only expected to give good results when the separation 
between two energy levels entered two-point function is large. In particular, it 
does not allow one to compute the correlations in the region where these are 
believed to be universal, namely for energy differences of the order of a level 
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spacing. 

The purpose of this paper is to compute expHcitly the two-point correlation 
function for the energy levels of the modular domain and the corresponding bil- 
liard. The calculations are based on a generalization of the Hardy-Littlewood 
method p3| which was originally developed to compute the distribution of prime 
numbers, and depend strongly on the number-theoretical properties of the mul- 
tiplicities of the periodic orbits of the modular group. 

This paper is organized as follows: In Sec. 2, we give a quick review of 
hyperbolic geometry and the Selberg trace formula and derive various basic 
relations between the spectrum (in particular its two-point function, which is 
the basic object of interest here) and the properties of the classical periodic 
orbits of the system. In Sec. 3, we describe a method originally due to Hardy 
and Littlewood to describe the structure of singularities of certain peculiar power 
series. This leads to the consideration of certain quantities which are evaluated 
exactly. These results are presented in detail in Sec. 4. This leads to a closed 
form for the two-point form factor of the spectrum of the modular domain which 
has (5-function singularities at all rational points. In Sec. 5, we demonstrate 
that after a suitable smoothing it becomes constant as it should for the Poisson 
distribution. This result is valid when the separation is fixed and energy tends 
to infinity. On a scale of hik/k (where k is the momentum) the two-point form 
factor has oscillations of number-theoretical origin. In Sec. 6 we generalize these 
formulas for the case of modular billiard with Dirichlct and Neumann boundary 
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conditions. In the conclusion (Sec. 7) we briefly repeat the main steps necessary 
for our derivation. In Appendices A-F the details of calculations are presented. 
Much standard material in hyperbolic geometry and number theory has been 
presented to make the article self-contained. 



2: Basic Identities 



In what follows we shall use the standard realization of the surface of constant 
negative curvature as the so-called Poincare plane, that is, the upper half-plane 
z — X + iy, where y > 0, endowed with the metric ds^ — y^'^{dx^ + dy^), (see 
e.g. ^). The hyperbolic distance between two points zi = xi + iyi and 

Z2 = X2+ iy2 is given by 

coshd(zi,Z2)^l + -^^1' (2.1) 
22/1^2 

The geodesies are then circles perpendicular to the real axis and the group of 
isometrics is the group of linear fractional transformations, that is 

, az + b ,„ „, 

z'^gz^^- (2.2) 
cz + a 

where a, 6, c and d are arbitrary real numbers, which can, without loss of 
generality, be chosen so as to satisfy the condition ad — be = 1. It is then easily 
verified that the composition of two such transformations gives another such 
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transformation according to the multiplication of the corresponding matrices 



One can now introduce finite surfaces through a device which has an analogue in 
the Euclidean case: There, one can construct a torus by identifying all points on 
the plane which differ from each other by a translation belonging to a discrete 
subgroup of the Euclidean group (which is the isometry group of the Euclidean 
plane). In the case of CNCS one proceeds similarly: One takes a given dis- 
crete subgroup G of SL{2, R) and identifies all points which are connected by a 
transformation g belonging to G. 

Under certain assumptions, this procedure leads to a CNCS on which free 
motion is ergodic. This motion has very strong chaotic properties |Q, ^, |^ and 
is therefore a natural object of study. 

To visualize this construction, it is convenient to introduce the notion of the 
fundamental domain of a given discrete group G. This is defined as a region in 
the upper half-plane such that 

1. No two points inside the fundamental domain arc connected by a trans- 
formation g belonging to G. 

2. For any point z' outside of the fundamental domain, there is a point z 
inside the fundamental domain such that there is a g in G with 




(2.3) 



z = 



gz. 



(2.4) 
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Eq. (2.4) leads to the identification of certain points on the boundary of the 
fundamental domain. Gluing these together yields a compact surface of constant 
negative curvature. The trajectory on the whole upper half-plane is given by 
a half-circle perpendicular to the real axis. After the identification has been 
carried out, this geodesic can be reduced to a curve lying entirely within the 
fundamental domain and consisting of segments of geodesies. 

Periodic orbits of this geodesic flow are in one-to-one correspondence with 
conjugacy classes of elements of G (see e.g. |^). If Af is any matrix belonging 
to G, then the (hyperbolic) length of the corresponding periodic orbit is given 
by the relation 

2cosh//2 = |TrM| (2.5) 

The natural "quantization" of such systems consists in the investigation of 
the spectrum of the invariant Laplace-Beltrami operator 

on the space of functions obeying the periodic boundary conditions 

*„(5^) = ^n{z) (2.7) 

for any element g of G. (Note that our definition of En differs by a factor 
of one half from the one commonly used in the literature.) Due to their pe- 
culiar mathematical structure, there exists for these models an exact relation 
between "quantum" eigenvalues and the periodic orbits of the corresponding 
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classical motion. This relation — the celebrated Selberg trace formula Q — 
can be stated in the following way: Let h(r) be an arbitrary even function with 
appropriate smoothness properties and g{u) its Fourier transform, given by 



Let the index p run over all classical periodic orbits, let Ip denote the length 
of the corresponding orbit and Lp be the length of the primitive periodic orbit. 
Let A be the area of the fundamental domain under consideration. One then 
has the identity: 



where r„ is equal to \j2En — 1/4, with i?„ being the eigenvalues of the Laplace- 
Beltrami operator with periodic boundary conditions. 

The first terms correspond to the smooth part of the level density and the 
second gives the contribution from periodic orbits. In mathematical language, 
it equals the sum over all conjugacy classes of matrices with trace larger than 
two. The "corner and horn" terms refer to contributions from group matrices 
with trace less than and equal to two respectively, whenever such elements exist. 
Their explicit form can be found in [^,0. 

In the following, we shall apply the Selberg trace formula to 




(2.8) 




OC 



+ "corner and horn" terms. 



(2.9) 



Mr)=<5(i?--(r2 + l/4)). 



(2.10) 
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This formally causes some problems, as the delta functions are not smooth 
enough for the series to converge. We will here proceed in a largely formal 
manner and will discuss the regularization later. 

In this way one obtains for the eigenvalue density d{E) — J2^=i ^ ^n) 

d{E) = {d{E)) + dosc{E) + d{E), (2.11) 

where the first term is given A/ (2^) 

do.c(i?) = E-TEy^n^- (2.12) 

ttA; 2smh/i,f^/2 

p n—l ^ 

and d{E) includes all other terms which enters the exact Selberg trace formula. 
Here k is the momentum, defined by E — fc^/2 + 1/8. 

Without the last term the Selberg trace formula agrees with the Gutzwiller 
trace formula. In principle, one could use a suitable regularization of the delta 
functions and take the appropriate limits at the end of the calculations. 

So far, the above formulae are valid for an arbitrary group. Now let us con- 
sider the simplest example of an arithmetic group, namely the modular group. 
This group is defined as the group of all 2 x 2 matrices with integer entries 
and unit determinant modulo the subgroup {1,-1}. It is well-known (see e.g. 
[pS)) that this group is generated by two of its elements: The translation T, 
which maps z to z + 1 and the inversion S, which maps z to — 1/z, with their 
corresponding matrices: 

-ill) -(-li)^ 
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The fundamental domain of the modular group has the form shown in Fig. 1 
and the area 7r/3. The arrows indicate the lines identified under the action of 



the generators T and S and the periodic boundary conditions of Eq. (2.6) mean 
that the function '^{x^y) takes the same values on corresponding lines. For 
this problem there is an evident reflection symmetry x —> —x, which leads to 
a splitting of the eigenfunctions in two classes, namely odd and even. These 
two classes correspond to the problem of finding the eigenvalues of the Laplace- 
Beltrami operator with Neumann and Dirichlet conditions on the boundary of 
the fundamental domain shown in Fig. 1. This problem is called the modular 
billiard (or Artin billiard |^) and we shall consider it in Sec. 6. 

We have said that the modular group is the group of all 2 x 2 matrices with 
integer entries and unit determinant. The periodic orbits on the modular domain 
then correspond in a unique way to the conjugacy classes of those elements of the 
modular group which have trace larger than two (hyperbolic elements) . Further 
one has the general relation 

|TrM| = 2cosh;p/2. (2.14) 

which connects the length of the periodic orbit Ip with the trace of a representa- 
tive matrix of the conjugacy class. But as all entries of the matrix are integers, 
the trace is also an integer. Here the arithmetical nature of the group comes 
into play. 

Therefore one sees that for all periodic orbits with a length less than L, the 
number of possible different lengths is the number of different integers N less 



than 2coshL/2, or asymptotically as L oo, N goes as e^/^. Now it is well 
known that for any group the number of periodic orbits of length less than L 
grows as § 0, H 

N{lp <L)^— (2.15) 

up to exponentially smaller terms. These two estimates show that in the case 
of the modular group (as well as other arithmetic groups) periodic orbits are 
degenerate, i.e., there are many periodic orbits with exactly the same length. 

Let g{l) be the number of periodic orbits of length /. The above estimates 
mean that asymptotically 



KL 

El = ^"-"^ (2-16) 

KL 

where the summation extends over different lengths of periodic orbits, counted 

without taking multiplicity into account. If we now define the mean multiplicity 

{g{l)) as follows 

, Number of periodic orbits with /< L < Z + A/ 

(q(1)) — - - (2 17) 

^' Number of different lengths with I < Ip < I + AV ^ ' ' 

one concludes that 

(5(0) = 2^. (2.18) 
For other arithmetic groups one obtains the same asymptotic behaviour but with 



a numerical prefactor which depends on the group |16, ^ This extraordi- 
nary degeneracy of the lengths of periodic orbits has been discussed before and 
is at the root of the remarkable structure to be found in these systems. 
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We now proceed to reexpress Eq. ( ^.11 ) in terms of a sum over all conjugacy 



classes of hyperbolic elements of the modular group. Denote by n the trace of 
a given conjugacy class and by g{n) the number of distinct conjugacy classes 
corresponding to trace n. Taking into account the fact that 7i goes as e'^/^ as 
n — )■ oo one concludes that 

2p^/2 n 
- -— - (2.19) 
L inn 

Therefore the mean multiplicity of periodic orbits having trace n grows asymp- 
totically as the number of primes less than n. While such a fact is suggestive 
of a deeper connection, the authors are as yet unable to state anything further. 
The multiplicity g{n) can also be identified with the proper class number of 
quadratic forms ([p5|, p6|), but we shall not require this representation in the 
following. 

Applying the Selberg trace formula to d{E), we split the contributions to 
dosc{E) into two parts. In the first we collect all periodic orbits whose matrices 
have traces less than a certain value ng 3> 1 and all others are putted in the 
second part. The first one we simply add to d{E) which from now will contain 
all "non-interesting" terms which are explicitly known and could be calculated 
without principal difficulties. Neglecting the difference between 2 cosh i and 

and further disregard the existence of multiple traversals (since they are 
exponentially few in number compared to the primitive orbits) we concentrate 
on 
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2 In 71 

dosciE) = — V 5(n)— cos(2fclnn). (2.20) 

TTK ^ — ' n 

n—UQ 

(All other terms are putted to d{E).) 

We have seen above that the coefficient entering in Eq. (2.20) is on the 
average of order one. Defining 

a{n) ^ g{n)\nnln, (2.21) 

we finally obtain that 

2 °° 

dosciE) = — V a(n)cos(2fclnn), (2.22) 

TTK 

n—riQ 

where {a{n)) is one. 

Eq. (2.22) as it is written diverges when the sum is performed over all values 
of n. This is the well-known divergence of the summation over long periodic 
orbits, which is inherent in all semiclassical formulae. Mathematically, one treats 
such problems by using a suitable function h(r) in the Selberg trace formula. 
We here use a different method of regularization based on the subtraction of the 
main term in the density of periodic orbits (see e.g., |l^). In Eq. (2.22), it 

corresponds to the substitution of a{n) by a{n) — a{n) — 1, i.e., the subtraction 
of the mean value of Q;(?i)|5 Of course, one has to add to d{E) the rest 

cos(2A:lnn) = 5RC(2iA:) - ^ cos(2fclnn), 

n—no n—1 

where (^(s) is the Riemann zeta function whose analytical continuation can be 
done easily by using the well-known functional relation . 

^This subtraction is not necessary for billiard problems with Dirichlct boundary conditions 
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Finally the level density of energy levels for the modular group can be written 
in the following form 

diE) = {d{E))+dosc{E), (2.23) 

where 

2 °° 

dosc{E) = — Q!(n) cos(2/c Inn), 

n=no 

and for simplicity we redefine {d{E)): 

{d{E)} = ^+d{E). 

Usually the mean density of levels includes only the Weil term (plus corrections). 
Sometimes it is convenient to put to it other terms as well. Here we shall consider 
all explicitly known terms and the convergent contributions from repetitions of 
primitive periodic orbits as part of {d{E)) . They are unessential for our purposes 
but can be important for numerical calculations. 

Having the expression for the level density, one can formally compute the 
n-point correlation functions: 

Rn{ei, . . . , e„) = {d{E + e^)d{E + e^) ■ ■ ■ d{E + e„)) . (2.24) 

The important point to note here is the energy smoothing denoted by the brack- 
ets: 

{f{E)) = J f{E')a{E - E')dE' (2.25) 

and / a{x)dx = 1. Here a{x) is a function which is peaked near x = and has 
a width AE. It is usually assumed that {d)~^ <C AE <C E, where (d) is the 
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mean level density. The standard choice is the Gaussian: 



<j{x) — exp 



but other smoothing fuctions are also possible. They should be choosen in such 
a way that if u is a constant (exp(ifcu)) —^ sufficiently fast as /c — > oo. This 
kind of averaging procedures is inevitable for the statistical analysis of a system 
in which there are no random parameters. 

In this paper we concentrate on the two-point function 

i?2(ei, £2) = {d{E + ei)d{E + e^)) , (2.26) 

where for the modular group we use the following expressions: 

d{E) = {d{E)) + dosc{E), (2.27) 



where as i? oo {d{E)) -^1/6 and dosc{E) is given by Eq. ( |2.22[ ). One then 
finds 

i?2(ei,e2) = {{diE)f) +i?2(ei,e2), (2.28) 

where 

i?2(ei,e2) = V a(ni)a(n2)(e2^('^^''^"-+^-^'""^) 

_^g2j(fcilnni-fc2ln«2) + ^.C.) (2.29) 



and ki is ■y/2(i? + e;) which goes as fc + e^/fc as — > oo. 

Due to the energy average, the first term will be washed out and the second 
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one will give contributions only when 712 = rii + r with r <C rii, 712- Finally 
Mei,e2) = ^ Oi{nMn + r) exp f 2z— - 2ie— j 



n=no r= — 00 



c. c. 



(2.30) 

where e = ei — e2. Note that the two-point function depends on the difference 
e = £2 — ei as it should be. 

Let us now assume that the following mean value exists 

1 ^ 

l{r) = ^lim^ XI '5(n)a(ri + r). (2.31) 

Then the dominant contribution to the two-point correlation function will be 

given by 

^2(e) = ;^3?y ^ 7(r-)e^'''^/"exp(^-2ze— j, (2.32) 

where we have used a continuum approximation for n, since only large values 

of n make a significant contribution. If we now define f{x) as follows: 

00 

m = E -^w^'"' (2-33) 

r= — 00 

we can finally express the two-point function as 

i?2(e) = — / e'="/V(2fce-'="/2)coseMdu, (2.34) 

where tq = 21nno/fc. 

Another frequently used quantity is its Fourier transform K{t), (the two- 
point form factor) 

1 foo 1 1 
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where w is 2ke-''*/^ and K{t) = when t < tq. Therefore all the non-trivial 
information is contained in the functions 7(r) and f{x). In the following section, 
we will outline a general method for evaluating them. The simplest approxima- 
tion, known as the diagonal approximation [ p2[ , would be to assume that the 
a{n) are essentially uncorrelated, that is, that 7(r) is zero for r 7^ 0. This gives 
for f{x) a constant value, which leads to an exponential growth of K{t) as is 
clear from Eq. (2.35) On the other hand, from a general consideration 

[ p2| (see Section 5) one sees that K(t) must saturate to a constant value for 
t ^ 00 if it was originally obtained from a discrete spectrum. The exact expres- 
sion for f{x) found in the next Section will give a resolution of this discrepancy, 
as will be seen in great detail in the final Sections. 

As a final remark, we note that Eq. (2.35) indicates that the modular domain 
in fact behaves much as an ordinary integrable system, in spite of its chaotic 
classical behaviour. Indeed, as we shall see, K{t) reaches the constant value of 
the order of one at a time t* « Ink/k which goes to zero as fc — > cx3. This means 
that in the region where the two-point correlation function is expected to be 
universal, namely for distances of the order of a mean level spacing, and hence 
for t not excessively small, the eigenvalues of the modular domain do not show 
correlation. On the other hand, there is structure present at small times, due 
to the cumulative effect of short periodic orbits. 



3: Two-point Correlation Function of Multiplicities 
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In this section we shall give a way of evaluating the quantities 7(r) and f{x) 
defined in the previous section. The first remark concerning the correlations 
7(r) is that they have undamped oscillations related to their number-theoretical 
nature. This recalls to some extent the correlation between prime numbers 
which show similar oscillations, as shown by the Hardy-Littlewood conjecture 
[ p3| . We shall therefore follow a similar line of approach to evaluate 7(r). 

The key point is to find a suitable expression for the multiplicity of periodic 
orbits of the modular group. In the previous Section we introduced the normal- 
ized multiplicities a{n) equal to g{n) In /n and it was shown that {a{n)) is equal 
to one. In Fig. 2 we present the correlations of the multiplicities a{n), that is, 
we show 7(r). What is striking about that graph is that it shows correlations 
which do not decay as r —^ oo. This is in a sense similar to the situation pre- 
vailing for the correlations between primes, as shown by the Hardy-Littlewood 
conjecture. We shall therefore follow a somewhat similar path. We define 



a(q; r) — lim — a(mq + r) 

rn— 

oc 

= lim(gM) ^ a(mg-|-r)e-(""'+'')", (3.1) 

m— 

where the equality between both is guaranteed by a Tauberian theorem 32, 33|. 
The intuitive meaning of this quantity is the average value of a[n) when n only 
runs over numbers of the form mq + r for given q and r. Clearly, since (a(n)) 



17 



is one 

9-1 

Y.a{q-r) = q. (3.2) 

Now, if a{n) were a smooth function, we would expect all a{q\r) to be equal 
to each other and hence to one. This, as we shall see, is not at all the case. 
Rather, the dependence of ck(g; r) on its arguments is exceedingly complex and 
highly irregular. This will in fact be one of the principal technical difficulties of 
this subject. 

In Appendix A, we shall show the exact way to compute a{q]r). To state 
the result, we need a few definitions: We call A/g the set of 2 x 2 matrices with 
entries being integers modulo q and having determinant one modulo q. These 
matrices form a group under multiplication modulo q and is sometimes called 
the modulary group |Q . Additionally, we define to be the set of elements 
of Mq with trace equal to r modulo q. We generally denote the number of 
elements of a set M by \M\. The result of Appendix A can then be stated as 
follows: 

= (3.3) 

In Table 1 we present the calculated values of a{q\r) for q < 11. These results 
are in fact in very good agreement with direct numerical computations of a{q; r). 
The intuitive meaning of this apparently strange result is the following: g{n) 
is the number of conjugacy classes of modular matrices of trace n. To each 
modular matrix, one can associate an element of Mg in a unique way simply by 
taking the entries of the matrix modulo q. If n is equal to r modulo q, then all 
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these matrices will belong to Mq^r- If we therefore assume that the matrices of 
the modular group cover the set Mq in some sense uniformly, Eq. (3.3) appears 
reasonable. The argument presented in Appendix A makes these ideas more 
rigorous. 

Now, in order to compute the correlation 7(r) between the a(n) we proceed 
as follows: Define as in psj 

oo 

$(z) = (3.4) 

Since {a{n)) is one, the convergence radius of this series is equal to one. The 
importance of this function comes from the fact that 

J,(e-") = e™ / (e-"+''^) $ (e-"-'^) e"*''^ = V a{n)a{n + r)e-2™ 

•^0 2^ n=l 

(3.5) 

and the right-hand side, again by a Tauberian theorem is connected to 

the quantity 7(r) which wc wish to obtain. 

The essence of the Hardy-Littlewood approach is the investigation of 
the function <i>(z) as z = exp(— u + + 2'Kip/q) as m — > and e — > 0, where p 
and q are coprime integers. The main step is then to write n in the form mq + r 
with r lying between and q — I and prove that in the expression for $(z) in 
Eq. (3.4) the dominant contribution as u and e go to zero will be given by the 
mean value of a{mq + r) , that is, one may substitute it by a{q; r). We present a 
more detailed discussion of the validity of this assumption in Appendix B. Here 
it may be sufficient to say that the basic assumption involved is the one that 
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there are no ordinary short-range correlations involved. Rather, all correlations 
have the oscillatory long-range behaviour shown in Fig. 2. Accepting this, one 
has that as m — > and e — > 

^ ^ ^ r=Om=0 

g-1 



- ^(^''^^ (3.6) 



u — le 



where 



0{P, q) = q-' E ^) (^ A , (3.7) 

r=o \ q / 

Therefore the function $(z) has a pole singularity at all points on the unit 
circle which have a rational multiple of 27r as phase. Its explicit form is given in 
Eq. (3.6). The reason for the appearance of such poles is the irregular behaviour 
of a(n) and, connected with it, the fact that the a{q; r) depend in a highly non- 
trivial way on q and r instead of being independent of them, as would be the 
case for more "reasonable" sequences. 

The next step is to substitute Eq. (3.6) into Eq. (3.5) We therefore divide 
the unit circle in intervals Ip^q centered around exp{2'jTip / q) , where p and q run 
over all relatively prime numbers with p < q and q less than some prescribed 
upper bound Q which later goes to infinity. If one now divides the integral in 
this way, neglects all terms in each interval except the pole terms of Eq. (3.6) 
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and finally extends the integration over e to the whole line, one obtains: 

i- l/?b,9)l^expf^r) (3.8) 



2u ^ ' V Q 

(P,9) = l ^ 

as w — !■ 0. Here and in the following, (p, q) will denote the greatest common 
divisor of p and q. Finally, from the definition of Jr(e^") in Eq. (3.5) one obtains 
for 7(r) 

7(0= E l/3b,9)l'expf^ry (3.9) 
The sum is performed over all 5, all < p < (7 and coprime to q, and the term 



p = and (7 = 1 is omitted as we defined 7(r) through aiji) (see 2.30) whose 
mean value is zero. 

This is the two-point correlation function of the multiplicities of the periodic 
orbits for the modular group . All other quantities of interest can be obtained 
from it. In particular the function j(x) introduced in the previous Section is 
given by 

/(x) = 2^ E |/3(p,g)|2<5(a;-27rp/g), (3.10) 

(P,9) = l 

where the summation is over all p and q coprime, without the restriction p < q. 

The subtraction needed to make the trace formula converge (see Section 2) 
is equivalent to setting /(O) equal to zero which is equivalent to removing from 
the sum in eq. (3.10) those terms for which p = mq, where m is an integer. 
From now on we assume implicitly that these terms are removed from the sum 
and the "renormalized" f{x) satisfies /(O) = 0. As we saw above, the knowledge 
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of f{x) determines immediately the function K(t) and consequently its Fourier 
transform, which is the two-point function. 

This is in essence the main result of our paper. Combined with a fairly 
technical evaluation of /3(p, q) which is carried out in Appendices C and D, this 
gives a closed form for the Fourier transform of the two-point function which is 
presumably exact. Note that in these evaluations there are no (or merely triv- 
ial) approximations. The substantial problems arise from the approximations 
involved in the use of the Hardy-Littlewood method as well as in the various 
simplifications of the Selberg trace formula. 

4: Results 

The fundamental equation worked out at the end of the preceding section 
gives an exact form for the two-point function as well as its Fourier transform, 
but they are still somewhat unwieldy. To simpliiy them and cast them in a 
useful form will be the purpose of this Section. 

To this end we need to point out a basic number-theoretical property of the 
functions we have been discussing: they are all so-called multiplicative functions. 
One says that the function g{n) is multiplicative when it has the following 
property 

g{mn) = g{m)g{n) whenever {m,n) = 1. (4.1) 
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Multiplicative functions are therefore uniquely determined by their values on 
numbers of the form where p is a prime number. It turns out that a{q; r) 
is multiplicative in the argument q and the functions /3{p, q) as well. The first 
follows from the fact that Mq^r and Mq can be expressed as the number of 
solutions of a given set of congruences modulo q. These are therefore multi- 
plicative in g as a consequence of the following well-known fact known as the 
Chinese Remainder Theorem(see e.g. p5|): If qi and (72 are relatively prime, 
then to every solution of a congruence or set of congruences modulo qiq2 there 
corresponds uniquely a pair of solutions to the same congruences mod qi and 
q2 respectively and vice versa. The fact that more complicated expressions such 
as q) maintain this multiplicative property is a tedious but straightforward 
exercise. 

To simplify the expression for /(x), we will use the following identity valid 
for an arbitrary multiplicative function g{n): 



00 



00 




n 




(4.2) 



p prime 



which is known as Euler's identity. This leads to 



00 



7(r) 



^A(n)-1 



00 




(4.3) 



p 



k=0 



where Ar{q) is given by 




(4.4) 
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To give a closed expression for Ar {q) we still need one standard definition from 
number theory: we define the Legendrc symbol of a and q, where q is an odd 
prime as follows 



= 1 if there is an x ^ (mod q) such that a — (mod q), 
= if a = (mod q) (4.5) 
= —1 otherwise. 

The meaning of this number is perhaps best understood by saying that the 
number of distinct solutions of the equation x'^ = a (mod g) is 1 + (a/q). 
The properties of the Legendre symbol are stated in any standard reference 
on number theory (see e.g. |3^, Q). A fairly tedious evaluation of Ar{q) 
(performed in Appendix D) gives the following for q a prime power (which is all 
that is necessary, since Ar{q) is also multiplicative in q): First, let q be equal 
to p", where p is an odd prime. Then we have for n = I: 



Arip) 



(p2 - 1)2 



'^g A^^-4)((x + r)2-4 | ^ 



(4.6) 



For 71 > 2, we have, letting t be an arbitrary non-zero number modulo p 

Arip" 



2(1 -1/p) r = (modp") 

1 J -2/p r = tp"-^ (mod p") 



p2«(l-p-2)] e(n,p)(l - 1/p) r = ±4 (mod p") 

~e{n,p)/p r = ±4 + tp"^^^ (mod p") 

(4.7) 

where e{n,p) takes the value —1 if n is odd and p is of the form 4fc + 3 and is 
equal to 1 in all other cases. For p = 2, we must list down individual cases for 
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low powers and eventually state a general rule: 



A(2) 
A (4) 

Ar{16) 

Ar{32) 



1/9 r = (mod 2) 

-1/9 r=l (mod 2) 

1/18 r = (mod 4) 

-1/18 r = 2 (mod 4) 



= 



1/(9-16) 







1 r = (mod 16) 
-1 r = 8 (mod 16) 



and finally, for the general case n > 6 

Ar{2") = l/(9-22"-4) < 
All the terms not explicitly shown above are of course equal to zero 



2 r = (mod 2") 

-2 r = 2"-i (mod 2") 

1 r = ±(4 + 2"-2) (mod 2") 

-1 r = ±(4 + 2"-2 + 2"-i) (mod 2") 



(4.8) 
(4.9) 
(4.10) 
(4.11) 

(4.12) 



(4.13) 



5: Two-point Form Factor of Energy Levels 



We have shown in Section 2 that the Fourier transform of the two-point 
correlation function of the energy levels for the modular group can be written 
as follows: 

m = ^fM, (5.1) 

where w = 2fce~''l*'/^ and 

f{x) = 27T mP,q)\^S{x-27Tp/q), (5.2) 

(p,9) = l 
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and (3{p, q) can be expressed in terms of Kloosterman sums as in Eq. (C.7). 

The two-point correlation function of any system with a non-degenerate 
discrete spectrum should have the following asymptotic behaviour p^ : 

R2{ei,e2)^ {d{E))S{e2-ei) as £3 - ei ^ (5.3) 

This is a simple consequence of the fact that, in the absence of systematic level 
degeneracies in the sum 

i?2(ei, £2) = 5] 5(£; - E,,, + e^)5{E - + e^) (5.4) 

ni ,712 

only the terms with iJ„j — En.^ are important in the limit 62 — > ei. From this 
the asymptotic relation of Eq. (5.3) follows. Consequently the two-particle form 
factor K{t) should have the following large i-asymptotics: 

K{t) ^ {d{E)) as t c» (5.5) 

(For the modular domain {d{E)) 1/6 as £^ ^ 00.) For the Poisson distri- 
bution, the form factor K{t) takes its asymptotic value throughout the range 
of t of order one (there is a non- universal region near t = 0). In the case of 
the standard matrix ensembles (GOE, GUE and GSE) these form factors reach 
their asymptotic value when t ^ 1. 

Although the above limiting behaviour follows from quite general consider- 
ations, it is generally not possible to derive it from the semiclassical formulae. 
Indeed, this fact depends on interference between non-degenerate long periodic 
orbits and is therefore quite difficult to ascertain. In the case of the modular 
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domain, our knowledge is by now sufficiently detailed that we can attempt an 
explicit verification of this behaviour. Of course, since the Selberg trace formula 
is exact and the spectrum of the modular domain is discrete, there is every rea- 
son to expect that it will work. Nevertheless, such a calculation provides an 
extremely powerful consistency check on the various approximations made in 
order to arrive at the final result. This is what we shall do in this Section. 



The "exact" form factor as given in Eq. (2.34) is really a sum of delta func- 
tions at the points 

tp.q = 7 In — . (5.6) 

k np 



Eq. ( ^.34 ) can then be rewritten as 



TT^fc 

(P,9) = l 



p 



S{t-tp,g). (5.7) 



The limit t oo then corresponds to the behaviour of f{x) as x — s- 0. 

In Fig. 3a we present a plot of f{x) computed from the periodic orbits of 
the modular domain using a direct Fourier transformation of the computed 7(r) 



with an appropriate smoothing procedure 

N 

N 



fc=0 

where the smoothing parameter A has been taken equal to 3/A^^ and N — 1000. 
The big peaks corresponding to small values of p and q are well pronounced, 
but its behaviour between peaks is not clear. It can be visualize if instead of 
f{x) one computes its integral 

Gix) = r f{y)dy. (5.9) 
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The plot of G{x) is given in Fig. 3b. These pictures suggest that it is convenient 
to separate the function f{x) into two parts 

f{x)^f<}:^{x) + fN{x), (5.10) 

where the first function when < a: < 27r is given by Eq. (|]^) but with finite 
number of terms for which q < N 

/^-(x) = 2^ J2 mp,q)\^6{x-2TTp/q), (5.11) 

(p,</)=l 

0<p<q<N 

and the second function fpj (x) includes contributions from all terms with q > N 
fN{x)=2TT ^ \p{p,q)\^6(x-2TTp/q), (5.12) 

(P,9) = l 
q>N 

The calculation of the first function is straightforward and only the second one 
needs a special attention. 

Of course, one can simply calculate /^"^(a;) for sufficiently large N but we 
shall show that adding to it an approximate expression for /jv (x) leads to much 
better approximation even for small N. 

It is shown in Appendix C that \f3{p, q)\ — > as q — > oo. Therefore all peaks 
in /Ar(x) are small and it is resonable to compute its mean value. We shall 
proceed in the following way (see e.g. [§). Let us define the function 

9n{x) = / fN{y)dy. (5.13) 
Jo 

If we now find a smooth function Gn(x) such that 

\GN{x)-gN{x)\^GN{x), (5.14) 
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then we shall say that gN^x) and GNix) are of the same order. We may then 
define 

{fN[x)) = — — — . (5.15) 

Note that this method is equivalent to smoothing the form factor over a small 
energy interval as in Eq. ( 2.25| ). 



Formally: 

gNix) = 27r J2 l/'b''?)!'' (5-16) 

(P,9) = l 
p/q<x/{2-K) 
q>N 

where the summation is over all coprime integers such that 

0<p/q<-^, q> N 

When q is fixed, q) is a number-theoretical function of p. It is in fact quite 
erratic. Therefore we conjecture that it can to some degree of approximation 
be replaced by its mean value over all p, which we call (3{q): 

Piq) = {\l3{p,q)\) = ^^ ] — , (5-17) 

2^(p,g) = l ^ 

where the summation is performed over all p coprime to and less than q. But 
it is easy to see that 

^ |/3(p,g)P=^(<z)=n^o«'), (5.18) 

(p,g) = l i 

where the Ui are the prime factors of q and rii the power with which they occur. 
These values have already been computed in Appendix D. It is worth pointing 
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out that in the case of an odd prime one can reduce the intractable expression 



of Eq. (4.5) to the foUowing: 



Mp) - ^^^2 - (5-19) 



Further, since one has, by elementary number theory |35| p4| , the following 
facts about the Euler function 



^(i)= E i = nfi-;^)' (5-20) 



(p,g) = l uj\q 

it follows that f3{q) is known explicitly (see Appendix F). Let 

GNix) = 27: J2 ^(9)- (5.21) 

(P:9) = l 

0<p/q<x/{2Tr) 
q>N 

The calculation of this function is performed in Appendix F. Here for clarity we 
consider a simple prototype of such function. 
Let 

Jn{x)= (5.22) 

0<p/q<x 
q>N 



It differs from the exact function (5.20) in two things. First, the sum includes 



all p and not only coprime to q. Second, the function f3{q) is substituted by its 
asymptotic behaviour. Now 

q=N+l ^ 

where [y] is an integer part of y. 
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When Nx <C 1 and N oo one can replace the sum by an integral and one 
obtains 



dq [q] 

q 



2 ^ n2 

71—1 



—x\ (5.24) 



When Nx ^ 1, one has [xq] ^ xq and 



q=N+l ^ 

Therefore at small x <C 1/N jffix) behaves as 7r^a;^/12 but at x (a; < 2n) 



it grows as x/N. The exact function (5.20) has the similar behaviour but the 



computations are more complicated. The details are presented in Appendix F. 
When x<$:l/N 

GNix) = (5.26) 

Therefore one has 

TTT 

{fN{x)) = — {x^O). (5.27) 

When x ^ 1/N and (27r — x) 1/N the asymptotics of Gn{x) changes. In 
this region 

GNix)-^^x, (5.28) 



where 



p{p + 1) 
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and the product is taken over all primes including p = 2. Numerically C ~ .704. 

Combining these values one concludes that the function (/at (a;)) as < x < 
27r approximately has the shape as in Fig. 4 and continues periodically beyond 
this interval. 

In our approximation the function f{x) can be written in the simple form: 



Together with the explicit expression for fff^ix) in Eq. ( 5.10D it gives a quite 



accurate description of the two-point form factor for the modular domain. 

In Fig. 5 we presented the difference between the 'exact' function G{x) com- 
puted by taking into account all terms with q < 1000 and the sum of these 
two terms for different value of N. Note the difference in scales with respect to 
Fig. 3. Quite good agreement is observed even for N — 20. 

The knowledge of the function (/Ar(x)) is particular important in the region 
of small X because it gives the dominant contribution to the asymptotics of the 
two-point correlation formfactor. From Eqs. ( 5.26| ) and (5.27) it follows that as 
a: ^ 

fix) -> ^x and Gix) ^x^. (5.29) 

In Fig. 6 we plot the function G{x) computed as the sum of all terms up to 
N = 500 and N = 1000 in the double logarithmic scale which is the most 



sensitive to small x behaviour. The solid line is the asymptotics (5.28) 



The most important consequence is that in the limit fc — > oo and t fixed the 
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two-point formfactor tends to the constant value 
as it should be for the Poisson distribution. 



(5.30) 



6: Billiard problems 



We have mentioned that eigenfunctions of the Laplace-Beltrami operator for 
the modular group can be classified by the parity with respect to the inversion 
X —X. The odd (even) functions are eigenfunctions of the billiard problem 
with the Dirichlet (Neumann) conditions on the boundary of half of the mod- 
ular domain (see Fig.l). In this section we compute the two-point correlation 
functions for these problems separately. 

In the modular billiard problem group matrices are 2x2 matrices with 
integer entries but with the determinant equals both 1 and —1 (see e.g. §). 
Matrices with determinant — 1 describe geometrically an inversion with respect 
to a circle and they correspond to the following transformation: 

, _ az* +6 



z 



cz* + d' 



(6.1) 



where z* is the complex conjugate of z. 

The periodic orbits of the billiard problems can be identified with classes of 
conjugated matrices (both with determinants ±1) but their length is given by 
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Eq. ( p.4[ ) only if the matrix determinant equals 1. If it equals —1 the length of 
corresponding periodic orbits should be computed by 

2sinh(//2) = |TrM|. (6.2) 

Physically matrices with determinant +1 (—1) correspond to periodic orbits 
with even (odd) number of reflections from the billiard boundary. 

For billiard problems there exist an exact Selberg-type trace formula which 
expresses the density of states through periodic orbits ^ . It looks like the 
usual trace formula ( p. 11 ) but periodic orbits with odd number of bounces have 



an additional minus sign for Dirichlet boundary conditions. 
Performing the same steps as in Section 2 one obtains: 

dba{E) = {diE)) + d{E) + dosc{E), (6.3) 

where {d{E)) ~ A/2tt is a smooth part of the level density. A here is the area 
of the modular billiard equals a half of that of the modular domain: A = 7r/6, 

dosc{E) = —ry] (a^^Hn) + ea^-\n)) cos{2k\nn). (6.4) 

TT/C 

n— no 

a^^^(n) here are the normahzcd multiplicities of periodic orbits corresponding 
to classes of conjugated matrices with determinants +1 and —1 correspondingly 

aW(n) =25^(71) — , (6.5) 
n 

where g^^\n) is the number of periodic orbits with determinant ±1 and trace 
equals n. The factor 2 is introduced for the convenience, e = — 1 for the 
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Dirichlet boundary conditions and e = 1 for the Neumann ones. The function 
d{E) contains all other terms. 

The arithmetic nature of the modular billiard group leads to the conclusion 
that 

;aW(n)) = Jim aW(n) = l. (6.6) 



n=l 

Similarly to Appendix A one concludes that quantities 

Af-l 

N'-^io N 



o^^'(g, r) = lim — a*-^' (mg + r) (6.7) 



m=0 

are equal to 

where MgJ^^ is the number of matrices with integer entries modulo q whose 
determinant equals ±1 and whose trace equals r, \Mg\ is the total number 
matrices modulo q with determinant 1. (The total number of matrices with 
determinant —1 will be the same.) 

After merely rephrasing arguments of previous Sections we find that the 
two-point correlation form factor of the modular billiard can be written in the 
following form: 

^(*) = T^if^^^H^) + 2/(+-^(«^) + f^-Hw)), (6.9) 

477 W 

where w = 2kex.p{—kt/2) and 



/(^^■^^)(a;) = 27r ^ fi(-^\p,q)(3^^^\p,q)d{x - 2n^), (6.10) 
(P,g)=i ^ 
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q 



-r 



where 

/3(^)(p,g) = -^a(^)(g,r)expf 

The exphcit formulas for [3'^~\p, q) are given in Appendix C. 

To compute the average behaviour of K(t) we have to know the average 
values of f''~~\x) and f^^ \x) as x — > 0. As in the previous Section one 
should first find the mean values of the product of two over all values of p: 

= ^ E P^''\p.q)P'^''\p.q)- (6.11) 

"^^^^ p:{p,q) = l 

The later sum is connected to functions AQ^''^'^\q) defined in Appendix D. If 
q = uj"^uj2^ . . . oj^'' is the canonical representation of q as the product of different 
primes oj then 

/?('-^^)(g) = ^n4'"'''('^r)- (6.12) 
Using this value it is shown in Appendix F that as x — > 

(/(-)(.x)) = (/(++)(.x)) = ^x, (/(+-)(x)) = 0(x3/2). (6.13) 

Therefore as i ^ Infc/fc 

m ^ ^, (6.14) 

which coincides with the Poisson value for this quantity {K{t) = A/{2ttY and 
A = 7r/6). 

Note that the last of relations ( 6.12| ) means that in the universal limit eigen- 



values of different symmetry classes (odd - even with respect to the inversion) 
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are uncorrelated. This property is usually taken for granted on "general con- 
siderations". But the problem does not seem to be trivial because the same 
periodic orbits enter the trace formulas for both odd and even states. Only 
their phases are different. A priori it is unclear that the cross term will van- 
ish. To our knowledge the modular billiard is the only dynamical system where 
the absence of correlation between states of different symmetry can be checked 
analytically. (On this subject see also pO|.) 

At small values oi t ^ Infc/fc the billiard form factor has peaks similar 
to the ones already found in the case of the full modular group. Note that the 
largest peak with p/q = 1/2 is absent in the billiard with the Dirichlet boundary 
conditions because —1 = 1 mod 2. 

7: Concluding Remarks 

In this paper we have computed the two-point correlation function for the 
energy levels of the modular group and the modular billiard. From the point of 
view of classical motion these systems are ergodic with strong chaotic properties. 
But their arithmetical nature leads to a very large degeneracy of lengths of 
periodic orbits and, as a consequence, to the fact that their two-point correlation 
function tends to the Poisson value typical for the integrable systems and not for 
chaotic ones. At large scale the two-point correlation function has prominent 
number-theoretical oscillations. 
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To clarify the (tedious) derivation we very briefly repeat the essential steps. 

• The Selberg trace formula allows to write the density of states as a sum 
over classical periodic orbits (the function {d{E)) is explicitly known): 

2 °° 

d{E) = {d{E)) + -r X! ct{n)cos{2k\nn), (7.1) 



where a{n) = a{n) — 1, 



irk 

n=no 



Inn 

a{n) = 9{n) — 



and g{n) is the number of periodic orbits with trace equals n. 



The two-point correlation form factor can be expressed as follows: 



K{t) = -^JH (7.2) 



where w = 2kexp{—kt/2), 
and 



+00 



1 ^ 

7(r) = lim — >^ a(n)de(n + r) 

^ N^oo N ^ ^ I ^ I 



n=l 



is the two-point correlation function for the multiplicities of periodic or- 
bits. 

Using a probabilistic approach we show that 

m=0 ' 3' 
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where \Mq\ is the total number of matrices with entries which are taken as 
integers modulo q and \Mq^r\ is the number of such matrices having trace 
equal to r modulo q. 

• A generalization of the Hardy-Littlewood method permits to find an ex- 
plicit expression for 7(r): 

l{r)= |/3(p,9)l'exp(27ri^r), (7.4) 

(P,9)=l ^ 

where the summation is taken over all q and all p < g coprime to q and 

9: 



1 

/3(p, 9) = - X! "(^5 exp(27ri-r). (7.5) 

^ r— ^ 

Introducing the Kloosterman sums 

S{n,m,c) = exp ( -^{nd + md~^) \ , 



c;=o 

/?(f^, q) can be written as follows: 

where w are the prime divisors of q. 

• These formulae give the explicit expression for the two-point correlation 

form factor 



TT^k 



-/3(p,a) 

p 



2 

J(t-tp,,). (7.7) 



where 

2 A;g 
tpg = 7ln — • 
k np 
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In the limit k oo and t fixed, the dominant contribution comes from 
terms with p/q <C 1. Smoothing over such values we show that in this 
limit K{t) has the constant Poisson value: 

Here A = n/3 is the area of the fundamental region of the modular group. 
For small t (of order of \nk/k) K{t) has number-theoretical oscillations 
due to cumulative contributions of degenerate periodic orbits. For very 
small values of t (of order of 1 /k) the two point form factor has 5 function 
peaks connected with short periodic orbits. 

Analogous formulas can be obtained for the modular billiards. As a byprod- 
uct we proved that the energy levels of different symmetry are uncorrelated. It 
seems that practically all such formulae can be generalized to other arithmetical 
groups. Though the modular group is by no means a generic system, it is the 
first ergodic dynamical system for which it is possible to compute explicitly the 
distribution of the energy levels. 
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Appendix A 



Here we want to show how the expression for a{q; r) is derived. One notes 
first that the modular group is generated by two elements s and t defined as 
follows: 




(A.1) 



Since s"^ = = 1, any element of the modular group can be represented as a se- 
quence of alternating s and t"' where a = ±1. Further, such representations are 
in fact unique. From this follows that to each conjugacy class there corresponds 
uniquely a word beginning with s and ending with a , up to cyclic permuta- 
tions. That is to say, cyclically equivalent words generate identical conjugacy 
classes, but no other words do. To this end, it is essential, however, to enforce 
the condition that the word begin with an s and end with a t'^. In the following 
we will therefore only consider words generated by the matrices 

mi = st m2 = st~^. (A. 2) 

Let us now define a{q; r) more precisely: Up to a normalization factor of q, it is 
the probability that a conjugacy class belonging to a given value of n should have 
a trace equal to r modulo g, after averaging over n. Here instead of averaging 
over n, we shall average over all conjugacy classes which require k symbols to 
generate and eventually consider the average over many and sufficiently large 
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values oi k. It seems highly reasonable to claim that these two averages should 
be equivalent. This is in fact the major assumption in this Appendix. 

It is obvious that there are exactly 2*^ different words generated by k symbols. 
Each of those words denotes a different conjugacy class, up to cyclic invariance. 
For the overwhelming majority of words, this simply introduces a fixed factor 
of k which docs not matter in this discussion. Those words for which the de- 
generacy factor is different from k are exponentially rare and can be neglected. 
We can thus view the words (or conjugacy classes) as arising from the following 
Markov process: Starting from the identity, at each step multiply the matrix to 
the right randomly either by mi or by TO2 with probability 1/2. In this way all 
conjugacy classes are generated with equal probability. If one now projects this 
Markov process down onto the set Mq by simply considering the entries of each 
matrix as integers modulo q, then the Markov process described above projects 
to a Markov chain on a finite space, namely Mq. Under such circumstances, very 
general theorems (p3|) guarantee the existence of an approach to equilibrium. 
For this purpose, one needs to show that the process is ergodic, that is, that 
every element of Mq can indeed be reached by an appropriate combination of 
mi and TO2. This is seen as follows: Consider an arbitrary element of Mq as an 
element of the modular group. As such, it can be represented as a product of 
s and t'^, though not necessarily of mi and 7712. To obtain the latter, note that 
there is a number k such that 

{stf = 1 (mod q). (A.3) 
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In fact k can be chosen equal ioq—\. Inserting this representation of the identity 
either before or after the word, one can always bring it into the required form. 

As a final point, we need to show that the equilibrium attained on Mq is 
indeed the uniform distribution. This is seen by noting that mi and TO2 are 
always diff'erent and invertible, so that the Markov process always connects an 
arbitrary matrix a with two different matrices ai and 02 with probability 1/2. 
Thus the probability of finding matrix a after k steps satisfies the equation 

Pk{a) = ^ (Pfc_i(amj-') +Pfc-i(am^')) • (A.4) 

This has the uniform solution as an equilibrium (fc-independent) solution. Since 
in a finite space the uniqueness of equilibrium is guaranteed, one can show that 
uniform distribution is indeed approached. Prom this follows the desired claim: 
Indeed, since a{q; r)/q is identified with the probability that a word fall upon a 
matrix of trace r modulo q, this probability is clearly given by the ratio of the 
number of matrices in Mq having trace r to the total number. 

There is still a slight caveat, however: It is well-known that cyclic behaviour 
cannot be excluded on general grounds. Indeed, odd-even oscillations are in 
fact observed for q equal to two. Nevertheless, it is readily seen that the eff'ect 
of such oscillations dies out when one averages over k. Thus, in the absence of 
oscillations, it would be sufficient to take the average over one set of conjugacy 
classes defined by a sufficiently large value of k. In the presence of oscillations 
we must still average over different values of k. In fact, it could probably be 
shown that cyclic behaviour for cycles larger than two cannot exist. It certainly 
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has never been observed up to now in the specific cases we have looked at. One 
might add that the transition matrix of the above Markov chain has some very 
remarkable properties: Its eigenvalues are highly degenerate and appear to lie 
on very specific loci of the complex plane. A great deal of this unexpected 
structure can be traced back to the fact that this matrix is invariant under the 
action of the modulary group (the modular group taken modulo q) and is related 
to the regular representation of the latter which is always reducible. We do not 
go any further into those details, however, because they are unnecessary to our 
immediate purpose. 



Appendix B 



Here we give another derivation of the Hardy-Littlewood method, which has 
the advantage of highlighting the approximations involved. We need to evaluate 
the integral on the left-hand side of eq. ( |3.4[ ) which expresses the function 7(r). 
We therefore divide the unit circle in intervals Ip q centered around exjp(2T:ip/ q) , 
where p and q run over all relatively prime numbers with p < q and q less than 
some prescribed upper bound Q which later goes to infinity. If one now divides 
the integral in this way and denotes the interval Ip^q shifted so as to be centered 
around z = 1 by Jp,q, one obtains: 

7(r) = lim(2w) ^ ^ a(n)a(n')e"("+"')" (B.l) 

(p,ij) = l n,n' 
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f dip f 2'iTip / ^ ■/ ' \ i\ 

/ — exp I (n — n — r) — i{n — n — r)ip 1 . 



To simplify the expression, one rewrites the sums over n and n' as sums over 
m',r' and m",r" respectively, where n = m'q + r' and r' is between and q—1, 
and similarly for n'. In this case, the integral over Jp^q is strongly oscillatory 
for m' ^ m" and one finds: 

2m ''"^ 
7(r) = lini — a{m'q + r')a{m'q + r") (B.2) 

^ (p,g)=l m' r,r'=0 

exp (-(2m'g + r' + r")u) e^^'''-''"-'"^ 

At this stage we make a fairly tricky approximation: In essence, we assume 
that there are no other correlations between the a{n) than those implied by the 
dependence on q and r of its average over numbers which are equal to r modulo 
q. In that sense, we rewrite eq. (7.2) as: 

7(r) = 9"^ X] H o:iQ'',r')a{q;r")exp (^^{r' - r" - r] 

[p,q)=^l r,r'—0 



1 

2 

<i — i /<->_•„ \ 

= E ^ 

(p,9)=l 



^ a(g; r') exp ( - 

r'=0 ^ 



exp[-^r). (B.3) 



Perhaps some examples of sequences a{n) for which the method works, and 
others for which it fails may be helpful. For example, let us consider the sequence 
ap{n) which is equal to 1 when n = (mod 2p) and —1 when n = p (mod 2p) 
where p is a given prime, and zero otherwise. It is easy to check that in this 
case a{q; r) is zero unless g is a multiple of 2p and r a multiple of p. The 
complete calculation of 7(r) for this sequence by the formulae given in the text 
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gives complete agreement with the exact result. On the other hand, let us now 
randomize this sequence in the following way: The numbers a((2fc + l)p) and 
a{2kp) are randomly interchanged with probability 1/2. The resulting sequence 
is clearly still correlated. In particular, ^{p) is non-zero. Nevertheless it is seen 
that all a{q; r) of this sequence are zero. Thus the Hardy-Littlewood method 
clearly fails to take into account mere short-range correlations. It is seen that 
adding this sequence to another, say the a{n) we have been studying in the text, 
clearly modifies the correlations 7(r) but in no way affects the a{q;r), so that 
one really must assume that there is no short-range correlations of this type 
present. Can this assumption be sustained for the a{n) we have been studying? 

At first sight, this appears very strange, since we are dealing with correlations 
over quite short ranges. However, from Appendix A, we see that the natural 
variable in which to study the development of correlations is not the trace 
n but rather the number of symbols necessary to generate a given conjugacy 
class. This number, however, does not vary smoothly at all with trace and 
the overwhelming majority of conjugacy classes with nearby traces have quite 
different number of symbols that generate them. Thus it is allowable to consider 
them as decorrelated as the Hardy-Littlewood method implicitly does. 



Appendix C 
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In the following we develop some basic tools to compute expressions for \Mq\ 
and |Mg,r|- The fundamental identity we shall be using is the following easily 
verified fact: 

,-igexp(^)=<5.,o, (C.l) 

where x is an integer between and q — 1, or else an arbitrary integer taken 
modulo q. We therefore find, for /3(p, q): 

^ r=0 \ 1 J 

9-1 



\Mq^r\ is the number of matrices where a,b,c,d are taken modulo q 

such that ad — be = 1 and a + d = r modulo q. 
Therefore 

1 9-1 9-1 



/ 2TTip 



r=0 a6cd=0 

Q— 1 q—1 q—1 



(C.3) 



-— ^ ^ ^ expf [m(arf- 6c- 1) +m'(a + d-r) +pr] j . 



0^ 

^ I ^' r=0a6cd=0mm'=0 

Note that, in the following, the Kronecker deltas must always be interpreted as 
being 1 if the two indices are equal modulo q. We keep the evaluation of |Mg| 
for the end and perform the sums over m' and r first. This gives: 

9-1 9-1 

bcd=Om=0 

Now we use the fact that {p,q) = 1 from which follows that the equation 
md = —p (mod q) cannot hold unless {d,q) = {m,q) = 1. Under these cir- 
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^ q-l q-1 dm \ 

= TTjT X! ^^md+pSi^'^'i>\ {dp-mhc-m)\. (C.4) 



cumstances, on the other hand, the equation has a unique sohition, namely 
m = —pd^^ (mod q), where the inverse is to be taken modulo q. It is well- 
known that if (d, q) = 1, this inverse exists and is unique. From this follows, 
after performing all sums: 



Note that such a formula cannot be used to obtain, say, \Mq\ by setting p to 
zero, since the assumption that (p, g) = 1 is essential in deriving eq. (C.5). 

This formula gives an explicit expression of the Fourier transformation of 
number of matrices modulo q with determinant equals 1. For billiard problems 
we shall need also the same quantity but for matrices with determinant — 1. 
Generalizing the previous arguments we obtain: 




(C.5) 




(C.6) 



where e = ±1 corresponds to matrices with determinant ±1. 



If one introduces the so-called Kloosterman sums [|6l 




(C.7) 



then /J^*^^ (p, q) can be written as 




S'(p,ep, q) 



(C.8) 



It is easy to see that /3*^^)(p, g) ^ as g ^ oo. More specifically. It has been 



shown that : 



|S'(n, m, c)| < d{c)^{m,n,c 




(C.9) 
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where d{c) is the number of divisors of c and (m, n, c) is the largest common 
divisor of m, n, q. Such estimates yield sharp bounds for the p{p, q). 

To finish, we need the value of \Mq\. This is found as follows: First, we can 
limit ourselves to the case of q = s", where s is a prime number, because \Mq\ 
is a multiplicative function of q by the Chinese Remainder Theorem. Consider 
first the case q equal to s. In this case, we are dealing with matrices over a field. 
The number of singular matrices is easily seen to be s'^ + — s. Indeed, there 
are — 1 different ways of choosing the first row to be a non-zero vector, and 
to each of those correspond s different vectors for the second row, which must 
be taken proportional to the first. If the first row is zero, on the other hand, 
any choice for the second row will yield a singular matrix. Thus the number of 
regular matrices is found to be (s — l)s(s^ — 1). These are distributed uniformly 
over s — 1 different non-zero values of the determinant, so that 

\M,\ = s{s'' -l) = s^{l-s-^). (C.IO) 

Now consider the case q = s^. Any numbers a, b, c and d satisfying ad — bc= 1 
(mod s") will also satisfy the same equation modulo s"~^. Let us then define 

a = ai+as''-^ (mods"), (C.ll) 

where ai is o taken modulo s"~^ and a is an appropriate number taken modulo 
s. Define similar numbers for b, c and d. One then has 

ad- bc = aidi - biCi + s'^~^{aiS + dia - bi-f - Ci/3) (mods"). (C.12) 



49 



Hence, to obtain a solution modulo s", we need an arbitrary solution modulo 
s"~^ and a set of numbers a, /3, 7 and S satisfying 

aiS + dia — bi'y — ci(3 = (mods). (C.13) 

But it is impossible that ai, 61, ci and di should all be zero. Therefore, let 
us assume, say, that ai 7^ 0. Then all possible sets of values for a, [3 and 7 
yield exactly one value for 5. Therefore to each solution modulo s"~^ there 
correspond exactly solutions modulo s", from which follows, for q = s" 

|M,| = |M,/,| = - (C.14) 

From this follows the general relation (see e.g. ||2l|, 

\M,\^q^\{{l-p-^). (C.15) 

p\q 

Appendix D 



In this Appendix, we give an expression for Ar{q) which is defined by the 
following relation 

Ar{q)= l/3b,9)l'exp('^r] . (D.l) 

Using the above expression for /3(p, q) one obtains: 

Mq) = T^^riq), (D.2) 



50 



where 



«r(5)= E E E expf^(rf + d-i-^-^-^-r-)). (D.3) 



p:(p,g) = l (d,g)=l (5,g) = l ^ ^ 

Now we can limit ourselves to evaluating ar{q) for q equal to s" and s a prime, 
since 0^(9) is in fact a multiplicative function of q. Let us first consider the case 
where q is equal to s and s is odd. In this case we have: 

ar{s) = gggexp(^(d + d-i-5-5-i-r) 

p=l d=l 5=1 ^ 



^^'''P{d + d-' -d-S-' -r)] (D.4) 



s 



5 — 1 S — 1 S— 1 

EEE-p( 

p=0 d=l 5=1 ^ 



Now one notices that the sum appearing in the last equation of Eq. (D.4) de- 
notes, up to a factor of s, the number of solutions of the equation 

d + d,-^ -S-6-'^ =r (mods). (D.5) 

One sees further that x can be represented in the form d + d~^ modulo s if and 
only if a;^ — 4 can be represented as a square. This follows from the fact that 
the equation 

d + d-'^ =x (mod s) (D.6) 

has the solution 

d=^(^x± - 4) (mod s). (D.7) 

Here again 1/2 must be taken modulo s, but since s is an odd prime this is 
always possible. Further, we see that the number of different representations of 
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X as d + d ^ is 1 + ((a;^ — 4)/g) , where we are using the Legendre symbol already 
defined in the text. Prom this we obtain the following representation of ar{s): 



_g/ (x'-4)((x + r)^-4) \ ^ 



a;=0 



where we have made use of the following standard identity 



x=0 



which we prove for completeness in Appendix E which gives some elementary 
properties of the Legendre symbol. For all practical purposes, this representa- 
tion is sufficiently explicit, and no more specific one could be found for general 
values of r. For r = and r = ±4, however, closed expressions can indeed be 
found: We shall only require the case r equal to zero. In this case the remaining 
sum involving Legendre symbols becomes s — 2, since these are identically equal 
to one except for x = ±2. From this we finally obtain for s an odd prime 



Ms) = ^r^^- (D-io) 



Now let us consider the case of q equal to s" In this case, we reexpress 
eq. (D.2) in the following way: Let Br{q) be the number of roots of the equation 

d + d-^ -6-6-'^ =r (mod 9). (D.ll) 
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In other words; 



q-l 




r) . (D.12) 



p=0 (d,9) = l {S,q)=l 



ar(s") was defined as a sum over p coprime to s". It is evident that it is 
equivalent to the sum over all p minus sum over p divisible on s. From this 
follows that we can express ar{q) in terms of Br{q): 



The factor in the second term comes from the fact that the summation over 
d and S in ar(s") runs from to s" — 1, whereas in Br{s"~^) it only runs from 
to s"-i - 1. 

We must therefore evaluate the number of roots of eq. (D.ll). There is an 
important simplification, however: Since we must eventually evaluate the differ- 
ence between the number of roots of eq. (D.ll) for s" and s times the number 
of those roots for s"~^, we can neglect all roots modulo s"~^ which generate 
exactly s roots modulo s", since these then cancel exactly. In the following we 
shall always associate to a root modulo s" the corresponding number modulo 
s"~^, which is also a root of the equation. 

This means that we can immediately discard all roots in which we do not 
have both d = ±1 (mod s) and S = ±1 (mod s). Indeed, let either of those 
two be different from 1 modulo s; for definiteness let it be d. In this case we 



o^(s") = 




(D.13) 
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have 

rf = di+as"-^ S = 5i+ps"--\ (D.14) 

where di is a number modulo s"^^, a is a number modulo s and similarly for di 
and p. In fact, we shall systematically use suffixes and Greek letters according 
to this convention. Putting this into eq. (D.ll) and developing the terms in 
s""^ as first order infinitesimals (which one may do since they have the same 
property modulo s" of being different from zero but vanishing in any power 
higher than the first) one obtains 

a(l - dj-^) - - 5f ^) = (mods). (D.15) 

From this follows that any value of /? determines a uniquely and that /3 can be 
chosen arbitrarily. Thus there are exactly s solutions modulo s" to every such 
solution modulo s"~^, so that only solutions with both d = ±1 (mod s) and 
5 = ±1 (mod s) need be considered. 

In this case we can express d and 6 in the following way 

n— 1 n— 1 

d = ±l + ^afcs'= 5 = ±l + ^/?feS^ (D.16) 
fe=i fe=i 

Again, we can put these expressions in the equation to be solved treating the 
expression as a formal power series which is cut off at order n. Two cases appear: 
Either all ak and f3k which would appear in terms higher than linear modulo 
s"~^ vanish or they do not. In the latter case, a straightforward generalization of 
the above argument shows that every solution modulo s"~^ generates s solutions 
modulo s", so that these can again be discarded. Finally, there remains the case 
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in which no terms but the linear ones are non-zero. In this case all these terms 
cancel, so that we are led to the conditions r = (mod s"~^) or r = ±4 
(mod s""^) for ar-(s") to be different from zero. 

We must now distinguish between the case where n is even or odd. Let us 
first consider the even case. We then define n to be 2k and consider first the 
case r equal to zero. This means we can set 

n—l n—1 

rf = ±l + ^a;s' S = ±l + ^/3is^, (D.17) 

l=k l=k 

where the signs of the leading ±1 must now be taken equal, so that r will be 
zero. In this case the solutions modulo s" are simply all possible choices of ai 
and Pi over the prescribed range, since they automatically cancel. That is, there 
are 2s" solutions, where the factor 2 comes from the possibility to choose both 
signs for the leading term. Similarly, one obtains for the same numbers taken 
modulo s"~^ that there are 2s"~^ solutions, from which follows 

ao(s") =2s"(l-s-i). (D.18) 

Similarly, for n even, it is easy to see that 

a±4(s") = s"(l - s"^), (D.19) 

since we have no possibility of choosing two signs. The rest of the argument is 
exactly as above, however. 

If we now consider the case where r is zero or ±4 modulo s"~^ but not 
modulo s", we find the following: It is not possible to have satisfy the equation 
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(C.IO) modulo s" with numbers of the form given by eq. (D.17). This imphes 
that for i ^ 0: 

atsn-^{s"-) = -2s"-s-\ (D.20) 

as well as 

a±4+t."-i(s") = -s"s"'. (D.21) 

Let us now consider n odd, that is, equal to 2fc + 1. In this case we must 
take into account all solutions of the form: 

n—l n—1 

d = ±l + ^a;s' 5 = ±l + ^As', (D.22) 

l=k l=k 

If we now compute from this expression, we find 

2k 
l=k 

Putting the expression in eq. (D.23) into the equation, we find that 

d + d-'-S-6-' = {al_, - (3l_,)s-\ (D.24) 

Therefore, the number of solutions of this equation is the number of zeroes of 
the expression a\ — 0^., which is 2s — 1 multiplied by the number of possible 
choices of the remaining parameters. We therefore have 

So(s^'=+^) - sBo(s2'=) = 2(2s - 1)52*= - 2s2fe+i. (D.25) 

Prom this one obtains again 

ao(s") = 2s"(l - s-i), (D.26) 
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so that the expression in the case of r equal to zero is unchanged. In the case 
of r equal to ±4, we must consider the number of times the expression + 
becomes zero. This either occurs 2s— 1 times or once only, depending on whether 
— 1 can be expressed as a square modulo s, which, as is well-known, depends 
on whether s is equal to 1 or —1 modulo 4. Thus, if s = 1 (mod 4) the above 
formulae are also recovered in the case of odd n, whereas in the case of s = — 1 
(mod 4) one obtains 



where e(n, s) is 1 either if n is even or s of the type 4m + 1 and —1 in the other 
case. 

It now remains to check the formulae in the only remaining case, which 
is when r is equal to zero or ±4 modulo s"~^ but not modulo s", still with 
n odd. Under these circumstances the entire reasoning shown above can be 
repeated word for word, except that at the point where we evaluate how often 
a1 — 0^. takes the value zero, we must now ask how many times it takes on a 
value different from zero. As we shall see, this last is independent of the value 
considered and is s — 1. This is verified as follows: Let N be the number of ways 
in which w can be expressed as the difference of two squares. One finds 



,2k _ g2fe+l 



(D.27) 



from which follows that generally speaking 



a±4(s")-e(n,6-)s"(l-s-^), 



(D.28) 



N 
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(D.29) 



where we have again used the identity (C.8) to simphfy the expression. From 
this we obtain 



(D.30) 



whereas for ±4 + ts^~^ we must ask how often the number aj, + 0f. takes on a 
fixed non-zero value. Again, depending on whether s = 1 (mod 4) or not, it 
takes the value s — 1 times or s + 1 times as is seen by looking at the sum 



x=0 



1 



(D.31) 



The remaining case s = 2 can be treated analogously. The results are presented 



in Eqs. (|4.8|)-([4.12|). 

The above-discussed values of Ar{q) correspond to matrices with determi- 
nant + 1. For billiard problems we need also matrices with determinant —1. 
The corresponding formulae for q) are presented in Appendix B. Instead 

of one function Ar (q) one has 3 functions 



p:{p,q) = l 



( 2T:ip 
1^' 



(D.32) 



Ai'^^\q) — Ar{q) and the two other functions can be computed in the same 
fashion as Ar{q). We omit the details of the calculations and present only final 
results. 

When g = s is a prime (s ^ 2) 



Ai^^'^^Hq) 



\M, 



x=0 



(x2-4ei)((a;-f r)2-4e2) 



- 1 



(D.33) 
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For later applications the value of Aq''^^ {s) is important. One obtains: 

4— )(s) = 4"^+) ifs = 4m + l 

= 4 ifs = 4TO + 3, (D.34) 

Ai+-\s) = ^^^^^-^^{sH,{l) ~ s ~ 1) (D.35) 

where 

ff,(fc) = Oif 

The explicit formula for -ffs(l) is not known but it is possible to show [Q 
that if (-1/s) = 1 

where a is a number for which (a/ s) = — 1. Therefore for any s 

\Hs{l)\ < Vs- 

For g = s" (n > 2) it follows that if (—1/s) = — 1 (i.e. there is no solution 
of the equation = ~1 mod s) than 

4— Hs") =4+-)(s") = 0. 

If (—1/s) = +1 so that there is a number x such that = — 1 mod s, then 

And finally, if r 7^ ±1 ± x, then a[+"^(s") = 0. 
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An important consequence of these relations is that for all s 

Appendix E 



In this Appendix we give some standard properties of the Legendre symbol. 
Since multiplication modulo g' is a group when g is a prime, one has the well 
known identity 

a;9-i = 1 (x^O) (E.l) 

From this follows that if a; is a square and is non-zero modulo q one has x^"^"^^/^ 
equal to one, whereas for x arbitrary it can only take on the values 1 and — 1. 
On the other hand, one finds that the squares of all positive numbers between 
1 and {q— l)/2 are in fact distinct, so that there are at least {q— l)/2 different 
squares. But a polynomial of degree n (with non-zero leading coefficient) cannot 
have more than n different roots. We therefore have the general result 

From this follows 




(E.3) 



as well as the result that —1 has a square root if and only if q is of the form 4k +1. 
Since we have exactly {q — l)/2 squares and the same number of non-squares 
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we find the simple result 



X 



x=0 

Further, one finds, for a ^ 0: 

9-1 



E ^ (E.5) 



.= 1 V 



(E.5) 



These are the only properties we have used in the text. Much more is known, 
however, and can be found in 134, pql. 



Appendix F 



In this Appendix we shall compute the function 

(P,9) = l 
0<p/ij<a;/(27r) 
q>N 

where (3{q) is defined as follows: For prime powers w" it has the value Ao(a;"), 



where Ao(p") is defined by Eqs. (4^) and ( |4.6| ) in the body of the text. (3{q) is 



then determined everywhere by the multiplicative property. Let 
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is a canonical representation of integer q into a product of primes uji then 

«"^?n^r— iKi— ^ 

where if oj 7^ 2 

1) = 1 - 2/w - n) = 2(1 - 1/w). 

For w = 2 

B(2, 1) = 1/4, B(2,2) = l/2, B(2,3) = 0, 

B(2,4) = 1, B(2,5)=0, B{2,n) = 2 if n > 6. (F.4) 

To take into account the inequalities p/g < x/2'k and q > N we find it convenient 
to use the following identity: 

itnL-t^ s \ 1 ifO<i.< 1 ' 

where e > 0. We therefore introduce the function 



(P,'Z) = 1 



One then obtains (ei, €2 > 0) 



Giv(x) = -i-,/ f/ fG(.,.') (F.7) 

The summation in 0(5, s') is taken over all q and all p coprime. The latter 
sum can be computed using the standard formula (see e.g. ]23[|) 

E /(p) = EE/(^'5H^)' (F.8) 

[p^q) — l k—1 5\q 
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where ^{6) is the Mobius function defined as a muhiphcative function which is 
zero on all numbers which are divisible by a square and satisfies 

m(1) = 1 flip) = -I 

for all primes p. This gives 

(p^q) — l k—1 6\q uj\q 

where C(s) is the Riemann zeta function and the product extends over all prime 
factors. 

Therefore 



(F.9) 



where tOi and Ui are as in Eq. (F.l). Let us define 



n— 1 



where k = 3 — s — s' . For u; 7^ 2 a direct computation gives 

/ ^ ,,,nk ^ ' ' 



n=2 



n—2 

1 / 1 LU^-^ 

^1-^-2- 



Therefore if 7^ 2 

1 uj^-^ - 1 
w'" — 1 



CfcH = l- — ~2 ^ ^ . (F.12) 
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Similarly 



We now rewrite Eq. (F.8) as a sum over all uji and rii. One can then perform the 



sums over the and obtains a factor Ck{ijJi) for each tOi. Finally one obtains 

(F.14) 

where the product is taken over all primes. 

This product converges when 1 < 5Rs < 2 — 5fts'. To obtain the fmiction G{x) 



one should compute the integral (F.6) along the line parallel to the imaginary 
axis real part of which lies in the above interval. Because TV ^ 1 it is possible to 
shift the contour of integration over s' right until it reaches the first singularity 
of G{s^ s'). It is easy to see that this singularity is a pole at s' = 2 — s coming 
from the product over all primes. Putting s' = 2 — s — e and assuming that 
e ^ one has 

But from Eqs. ( |F.11| ) and ( |F.12D it follows that 



CiH = l-a;-2 (F.16) 
for all UJ including uj — 2. Therefore 

(.1 ^ 



64 



where 



and as e ^ 



G{s, 2-s-e)^{^Y ^"'C(s)C(l + ^)Ks. (F.18) 



2tt I 



where 



because 



id 

It is weh known that the Riemann zeta function C,{s) has a pole at s = 1 with 
unit residue (see e.g. [^) and consequently as s' ^ 2 — s 

Integrating over s' one concludes that in the leading order of 1 /N 

For convergence it is necessary that 1 < e < 2. Rewriting {xNy as exp(— s ln(a;iV)) 
one easily concludes that if xN <^ 1 one can shift the contour of integration 
right up to the pole at s = 2 and 

G^(x) ^ ^ (F.21) 

because ^(2) = 7r^/6 and K2 = 1. If xN ^ 1 one can move the contour of 
integration only left up to the pole s = 1 coming from C(s) and 

Gn{x) ^K, (F.22) 
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where 

K^ = --)(! + -IT^—T,) = 11(1 - 
and the product is taken over aU primes. 

In Fig. 5 we present the plot oi g{x). On a double logarithmic scale, the 
behaviour for x ^ 1 is in evidence. The prefactor is also found to be correct, as 
shown by the theoretical line. As indicated in the figure captions, we then show 
g{x) in which the peaks with denominators less than N have been subtracted. 
In these one sees an approximately linear behaviour over a broad range of x, 
which is in complete agreement with the behaviour predicted in Sec. 5 and in 
this Appendix. 

For billiard problems one has to compute functions G{x) with [3'^'^^ ''^'^\q) 



defined in Eq. (3.11). Using the values of A\^^''''^\q) presented in Appendix D 



one concludes that 

if (— l/w) = 1 and consequently 



for such primes in particular (see Eq. (F.15)) 

C['-\lo) = 1- 

For primes with (— l/tj) = — 1 

B'^—\uj,l) = l-uj-'^, B^—\uj,n)^Q if n > 2, 
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It means that for all primes the value of c[ (uj) is the same as for c[^^^ (uj) 
and the asymptotics of the function G^^^^x) coincides with that of G'^++^(.t). 
In particular as x ^ 



G(—\x) 



X 
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and when x':^ 1 /N 



For the function •* {x) one obtains from Eqs. (|D.34| ) 
B^+-\u,l) = \{sH,{l)-s-l) 
and B^+-\uj, n) = for n > 2. Therefore 

G\;-\u:) = \{sH,{l)-s-l). 

UJ 

But we have mentioned that 



(F.23) 



(F.24) 



\li,{\)\<^s. 



It means that the product (F.13) will have a singularity only when s — 2.5 and 

as a; < 1 /A^ 

g'^^-\x)=0{x'-''). (F.25) 
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Figure captions 

Fig. 1 The fundamental domain of the modular group. T denotes the transla- 
tion z ^ z + 1, S IS the inversion z — > — 1/z and the arrows connect the 
boundaries identified under these transformations. The shaded region is 
the fundamental domain of the modular billiard. 

Fig. 2 Two-point correlation function of the multiplicities of the periodic or- 
bits of the modular group computed from the knowledge of whose corre- 
sponding matrices have traces up to 8000. 

Fig. 3 a) The Fourier transform of the two-point correlation function of the 
multiplicities of the periodic orbits for the modular domain. 

b) Its integral. 

Fig. 4 The schematic picture of the function (/at (a;)) for < a; < 2tt. 

Fig. 5 The difference between the exact G{x) and the approximate formula 



(5.28) for different values of N which are indicated near the curves. The 



middle line corresponds to iV = 50. 

Fig. 6 Behaviour of the integral of the two-point correlation formfactor at small 
X. The dotted line corresponds to the sum of all terms up to iV = 500 and 



the solid one to = 1000. The straight line indicates asymptotics (5.26). 
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